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Abstract 

In this paper, we compute the closed form expressions of elastody- 
namic Green functions for three different viscoelastic media with simple 
type of anisotropy. We follow Burridge aZ. [Proc. Royal Soc. of London. 
440(1910): (1993)] to express unknown Green function in terms of three 
scalar functions by using the spectral decomposition of the Christoffel 
tensor associated with the medium. The problem of computing Green 
function is, thus reduced to the resolution of three scalar wave equations 
satisfied by 0i, and subsequent equations with (fii as source terms. To 
describe viscosity effects, we choose an empirical power law model which 
becomes well known Voigt model for quadratic frequency losses. 
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1 Introduction 

Numerous applications in biomedical imaging [6, 14], seismology [2, 22], ex- 
ploration geophysics [29, 30], material sciences [4, 15] and engineering sciences 
[1, 17, 32] have fueled research and development in theory of elasticity. Elas- 
tic properties and attributes have gained interest in the recent decades as a 
diagnostic tool for non-invasive imaging [28, 37]. Their high correlation with 
the pathology and the underlying structure of soft tissues has inspired many 
investigations in biomedical imaging and led to many interesting mathematical 
problems [7, 10, 9, 11, 8, 16, 38, 39]. 

Biological materials are often assumed to be isotropic and inviscid with re- 
spect to elastic deformation. However, several recent studies indicate that many 
soft tissues exhibit anisotropic and viscoelastic behavior [27, 35, 38, 39, 33, 47]. 
Sinkus et al. have inferred in [38] that breast tumor tends to be anisotropic, 
while Weaver et al. [46] have provided an evidence that even non cancerous 
breast tissue is anisotropic. White matter in brain [33] and cortical bones [47] 
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also exhibit similar behavior. Moreover, it has been observed that the shear 
velocities parallel and orthogonal to the fiber direction in forearm [35] and bi- 
ceps [27] are different. This indicates that the skeletal muscles with directional 
structure are actually anisotropic. Thus, an assumption of isotropy can lead 
to erroneous forward-modelled wave synthetics, while an estimation of viscosity 
effects can be very useful in characterization and identification of anomaly [16]. 

A possible approach to handle viscosity effects on image reconstruction has 
been proposed in [18] using stationary phase theorem. It is shown that the ideal 
Green function (in an inviscid regime) can be approximated from the viscous 
one by solving an ordinary differential equation. Once the ideal Green function 
is known one can identify a possible anomaly using imaging algorithms such 
as time reversal, back-propagation, Kirchhoff migration or MUSIC [7, 12, 14, 
6]. One can also find the elastic moduli of the anomaly using the asymptotic 
formalism and reconstructing a certain polarization tensor in the far field [10, 
12, 15, 13]. 

The importance of Green function stems from its role as a tool for the numer- 
ical and asymptotic techniques in biomedical imaging. Many inverse problems 
involving the estimation and acquisition of elastic parameters become tractable 
once the associated Green function is computed [5, 7, 12, 18]. Several attempts 
have been made to compute Green functions in purely elastic and/or isotropic 
regime. (See e.g. [18, 16, 19, 22, 36, 43, 44, 45] and references therein). How- 
ever, it is not possible to give a closed form expression for general anisotropic 
Green functions without certain restrictions on the media. In this work, we pro- 
vide anisotropic viscoelastic Green function in closed form for three particular 
anisotropic media. 

The elastodynamic Green function in isotropic media is calculated by sep- 
arating wave modes using Helmholtz decomposition of the elastic wavefield 
[2, 18, 16]. Unfortunately, this simple approach docs not work in anisotropic 
media, where three different waves propagate with different phase velocities and 
polarization directions [22, 17, 23]. A polarization direction of quasi-longitudinal 
wave that differs from that of wave vector, impedes Helmholtz decomposition 
to completely separate wave modes [26] . 

The phase velocities and polarization vectors are the eigenvalues and eigen- 
vectors of the Christoffcl tensor F associated with the medium. So, the wavefield 
can always be decomposed using the spectral basis of F. Based on this observa- 
tion, Burridge et al. [19] proposed a new approach to calculate elastodynamic 
Green functions. Their approach consists of finding the eigenvalues and eigen- 
vectors of Christoffel tensor r(Va;) using the duality between algebraic and 
differential objects. Therefore it is possible to express the Green function G in 
terms of three scalar functions 4>i satisfying partial differential equations with 
constant coefficients. Thus the problem of computing G reduces to the reso- 
lution of three differential equations for (j)i and of three subsequent equations 
(which may or may not be differential equations) with (f>i as source terms. See 
[19] for more details. 

Finding the closed form expressions of the eigenvalues of Christoffel tensor 
r is usually not so trivial because its characteristic equation is a polynomial 
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of degree six in the components of its argument vector. However, with some 
restrictions on the material, roots of the characteristic equation can be given 
[36]. In this article, we consider three different media for which not only the 
cxplicite expressions of the eigenvalues of T are known [19, 44], but they are also 
quadratic homogeneous forms, in the components of the argument vector. As a 
consequence, equations satisfied by become scalar wave equations. Following 
Burridge et al. [19], we find the viscoelastic Green functions for each medium. 
It is important to note that the elastodynamic Green function in a purely elastic 
regime, for the media under consideration, are well known [44, 19]. Also, the 
expression of the Green function for viscoelastic isotropic medium, which is 
computed as a special case, matches the one provided in [18]. 

It has been shown in [20] that Voigt model is well adopted to describe the 
viscosity response of many soft tissues to low frequency excitations. In this 
work, we consider a more general model proposed by Szabo and Wu in [40], 
which describes an empirical power law behavior of many viscoelastic materials 
including human myocardium. This model is based on a time-domain statement 
of causality [41, 42] and reduces to Voigt model for the specific case of quadratic 
frequency losses. 

We provide some mathematical notions, theme and the outlines of the article 
in the next section. 

2 Mathematical Context and Paper Outlines 

2.1 Viscoelastic Wave Equation 

Consider an open subset ft of R"^, filled with a homogeneous anisotropic vis- 
coelastic material. Let 

u(x,t) : X K+ ^ 

be the displacement field at time t of the material particle at position x € 51 
and Va;u(x,t) be its gradient. 

Under the assumptions of linearity and small perturbations, we define the 
order two strain tensor by 

£ : (x, i) e 17 X R+ i (V,u + V,u^) (x, t), (2.1) 

where the superscript T indicates a transpose operation. 

Let Q e and Y e >Cs(M^) be the stifi'ness and viscosity tensors of 

the material respectively. Here £s(K'^) is the space of symmetric tensors of 
order four. These tensors are assumed to be positive definite, i.e. there exists a 
constant S > such that 

{Q-e)-i>s\e and {^■i)-i>s\e, ve 

where £s(M'^) denotes the space of symmetric tensors of order two. 
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The generalized Hooke's Law [40] for power law media states that the stress 
distribution 

CT : O X M+ ^ CsiM.^) 



produced by deformation e, satisfies: 

SL^Q-£+X- A[e] 
where ^ is a causal operator defined as 



(2.2) 



_(_l)7/2 



dt-y- 



7 IS an even integer, 



-(7 - ^ is an odd integer, 



(2.3) 



--r(7) sm(77r/2)-^*t^ 



7 is a non integer. 



Note that by convention, 

A[u]i = A[ui] and .4[e]ij = A[eij] 1 < «, J < 3. 

Here H(t) is the Heaviside function, F is the gamma function and *t represents 
convolution with respect to variable t. See [3, 21, 40, 41, 42] for comprehensive 
details and discussion on fractional attenuation models, causality and the loss 
operator A. 

The viscoelastic wave equation satisfied by the displacement field u(x, t) 
reads now 

where F(x, t) is the applied force and p is the density (supposed to be constant) 
of the material. 

Remark 2.1 For quadratic frequency losses, i.e, when 7 = 2, operator A re- 
duces to a first order time derivative. Therefore, power-law attenuation model 
turns out to be the Voigt model in this case. 

2.2 Spectral decomposition by Christoffel tensors 

Wc introduce now the Christoffel tensors F*^,!^^ : M''* £s(M'^) associated 
respectively with C and V defined by: 



rzj(n) = C^kUjUkUj, r^j(n) = ^ VkUjUkUj, 



Vne 



1 < t,j < 3. 



k,l=l 
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Remark that the viscoelastic wave equation can be rewritten in terms of Christof- 
fel tensors as : 

P^-F=r[V.]u + r[V.M[u]. (2.4) 

Note that r*^ and arc symmetric and positive definite as C and V are already 
symmetric positive definite. 

Let L'i be the eigenvalues and be the associated eigenvectors of T'^ for 
i = f , 2, 3. We define the quantities Mf and by 

= D^D^, and = (Ai-)-iD? (2.5) 

As r'^ is symmetric, the eigenvectors are orthogonal and the spectral 
decomposition of the Christoffel tensor F'^ can be given as: 

r = E^»'E.^ ^ith I = ^E^ (2.6) 

i=l 4=1 

where I e £s(M'') is the identity tensor. 

Similarly, consider F" the Christoffel tensor associated with V and define 
the quantities L^, D", and E^ such as 

3 3 
r" = ^LrEr with (2.7) 

4=1 1=1 

We assume that the tensors F'^ and F" have the same structure in the sense 
that the eigenvectors and are equal. (See Remark 3.3). In the sequel we 
use D instead of D'^ or D" and similar for E and M, by abuse of notation. 



2.3 Paper Outline 

The aim of this work is to compute the elastodynamic Green function G asso- 
ciated to viscoelastic wave equation (2.4). More precisely, G is the solution of 
the equation 

(rlV,]G(x,t) +r [V.]^[G](x,t)) - P^^=^ = Sm^)l, (2.8) 
The idea is to use the spectral decomposition of G of the form 

3 3 

G = $^E,(V,)<^,; = ^(D,; ® D,)A/-V., (2.9) 

i=l i=l 

where 4)i are three scalar functions satisfying 

(L?(V,)0. + mV,)A[^.,]) - = <5(t)5(x) (2.10) 

(See Appendix A for more details about this decomposition.) 
Therefore, to obtain an expression of G, we need to: 
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1- solve three partial differential equations (2.10) in (j>i 

2- subsequent equations 

V', = MrVi (2.11) 

3- and calculate second order derivatives of ipi to compute 

(Di (g) D,)7/'i 

In the following Section, we give simple examples of anisotropic media which 
satisfy some restrictive properties and assumptions (see Subsection 3.4) defining 
the limits of our approach. In Section 4, we derive the solutions 0i of equations 
(2.10). In Section 5, we give an exphcite resolution of tpi = M~^(j)i and (Dj (g) 
'Di)ipi. Finally, in the last section, we compute the Green function for three 
simple anisotropic media. 



3 Some Simple Anisotropic Viscoelastic Media 

In this section, we present three viscoelastic media with simple type of anisotropy. 
We also describe some important properties of the media and our basic assump- 
tions in this article. 

Definition 3.1 We will call a tensor c ~ (c„i„) G £s(R^) the Voigt represen- 
tation of an order four tensor C G £g(R'^) if 

Cmn — Cp^ij^p(^i^j'j = Cijkl 1 ^ i, J, fcj ^ ^ 3 

where 

p{i,i)=i, p{i,j)=p{j,i), p(2,3)=4, p(l,3) = 5, p(l,2) = 6. 

We will use c and v for the Voigt representations of stiffness tensor C and 
viscosity tensor V respectively. 

We will let tensors c and v to have a same structure. For each media, 
the expressions for T", L-(V^), and M^{\7^,) are provided [19, 44]. 

Throughout this section, will assume the value Cpg for c and Vpq for v 
where the subscripts p,q £ {1, 2, • • • , 6}. Moreover, we assume that the axes 
of material are identical with the Cartesian coordinate axes 61,62 and 63 and 
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3.1 Medium I 



The first medium for which we present a closed form elastodynamic Green func- 
tion is an orthorhombic medium with the tensors c and v of the form: 
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^M55 
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M66/ 



The Christoffcl tensor is given by 

/ Cuff! + ceedl + C55di 

r = ceed^ + csa^l + Ci^di 

\ Css^i^ + C44a| + C33i9| 

Its eigenvalues L'^{\7x) and the associated eigenvectors D'i{\7x) are: 
i^(V,) = cndf + ceedl + 0558^ 

L2C^x) = C66(9i + 02282 + C44di 
L^aiVx) = C55d'l + C44a| + C33dl 

= e, with = 1 Vi = 1, 2, 3 



3.2 Medium II 

The second medium which we consider is a transversely isotropic medium having 
symmetry axis along 63 and defined by the stiffness and the viscosity tensors c 
and V of the form: 
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M66/ 



with /i66 = (mii - Mi2)/2. Here 

/ ciidf + cqqOI + Ciidl (cii - cee) ^i^a 
= {cii - cqq) 8182 CQedl + ciidl + didl 

\ C449? + C449| + C33C)| 

The eigenvalues Lf(Vx) of T^iV x) in this case are 

LliVx) - Ciidl + C44a| + C339| 
i2(Vx) =Cii52+Ciia|+C44a| 

Ll{\/x) = ceea^^ + £66^2 + C44af 
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and the associated eigenvectors 'D'i{Vx) are: 




D2 = 

Thus Mf = 1, and M| ^ M§ = df 
3.3 Medium III 

Finally, wc will present the elastodynamic Green function for another trans- 
versely isotropic media with the axis of symmetry along 63 and having c and v 
of the form: 



/ Mil 


//II — 2/i66 


/111 — 2/144 
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/ill 


/ill — 2/144 
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The Christoffel tensor in this case is 



' ciid{ + ceed^ + Ciidi 
(cii - cqg) did2 
(cii - C44) didz 

Its eigenvalues Lliy^) are: 

iJ(V,) = 



and the eigenvectors Tiliy x) are: 
/ai\ 

= I 92 1 , d;5 



(cii - CQ(i)did2 
ceedf + cii9| + C445| 

(Cii - C44) ^293 



(cii - C44)9ia3 
(cii - C44) 92^3 
C449^ + C449I + Cii9| 



H Cii92^ - 


h cii9f 


= C11A3; 


1- C6692 - 


h C4493 




1- C449I - 


h C449| 


= C44A^ 


92 \ 




/ -9i93 




U3 — 


-9293 









(3.1) 



In this case, Mf = = df + 9| and M§ = {df + 9|)A^ 

3.4 Properties of the Media and Main Assumptions 

In all anisotropic media discussed above, it holds that 

• The Christoffel tensors T_'^ and F" have the same structure in the sense 
that 

D'^^D^, Vi = 1,2,3. 
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• The eigenvalues L°(Va;) are homogeneous quadratic forms in the compo- 
nents of the argument vector V^; i.e. 



92 

and therefore equations (2.10) are actuahy scalar wave equations. 
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• In all the concerning cases, the operator (V^r) is either constant or has 
a homogeneous quadratic form 
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In addition, we assume that 

• the eigenvalues of and satisfy 

L,''(V.) = AL,^(V.). 



• and the loss per wave length is small, i.e. 



ft « 1. 



Remark 3.2 The expression M| = (pl + c^DA^, wiZ/ 6e avoided in the con- 
struction of the Green function by using the expression 

G = 031 + El (V,) (01 E2 (V,) (02 - 03) 

for the elastodynamic Green function. 

Remark 3.3 In general, and are dependant on the parameters Cpq and 
Vpq . Consequently, Tf' and T_^' can not be diagonalized simultaneously. However, 
in certain restrictive cases where the polarization directions of different wave 
modes (i.e. quasi longitudinal (qP) and quasi shear waves (qSH and qSV)) are 
independent of the stiffness or viscosity parameters, it is possible to diagonalize 
both r'^ and simultaneously. Moreover, the assumption on the eigenvalues 
and LI, implies that for a given wave mode, the decay rate of its velocity in 
different directions is uniform, but for different wave modes (qP, qSH and qSV) 
these decay rates are different. 



4 Solution of the Model Wave Problem 

Let us now study the scalar wave problems (2.10). We consider a model problem 
and drop the subscript for brevity in this section as well as in the next section. 
Consider 

(L^[v,]0 + L^[vMm - = Km^)- (4.1) 
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Our assumptions on the media imply that L'^ and L" have the following 
form; 

ilVj=^a|^ and L"[V.] = /3L^[V.] = ^ /3af ^ 



1 ' ^^i 



Therefore, the model equation (4.1) can be rewritten as: 



-P0 = <5(t)<5(x), 



By a change of variables Xj = we obtain in function (/)(^) = (^(x) the 

following transformed equation : 



A, 



(4.2) 



where the constant a = 010203. 

Now, we apply A on both sides of the equation (4.2), and replace the result- 
ing expression for A A^cj) back into the equation (4.2). This yields: 



A^cj) + j3A 



92 



9*2 



Recall that /3 << 1 and the term in /32 is negligible. Therefore, it holds 



■13 A 



9*2 



9*2 a 



(4.3) 



Finally, taking temporal Fourier transform on both sides of (4.3), we obtain the 
corresponding Hclmholtz equation: 



+ ( 1 - pA{uj) ) $ = ( 1 - pA{io) ) ^5{0 



(4.4) 



where $(^,a;) and A{uj) are the Fourier transforms of 4'{(,,t) and the kernel of 
the convolution operator A respectively. Let 



k{uj) = ^^2 i^i- i3A{uj)j. 

Then the solution of the Hclmholtz equation (4.4) (see for instance [25, 34]) is 
expressed as 

g%/^K(w)T(x) 



<i>{^,Lj) = ^[l-(3Aiu;) 



4a7rT(x) 
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where 



/ tJC^ 



Using density normalized constants = — ^, we have 

'VP 



$(x,.)Hl-MH)^^-^. (4.5) 



where constant h = 6162^3 and 



/ 'y»2 rv*^ 

5 Solution of the Model Potential Problem 

In this section, we find the solution of equation (2.11). We once again proceed 
with a model problem. Once the solution is obtained, we will aim to calculate, 
its second order derivatives for the evaluation of D (g) Di/;. 

5.1 Solution of the Potential Problem 

Let V'(x, t), be the solution of equation (2.11) and 4'(x, w) be its Fourier trans- 
form with respect to variable t. Then \I'(x, w) satisfies, 

= ^(x,-) = (1 - m^)) "^ll'^^l^^ ■ (5-1) 

When M is constant, the solution of this equation is directly calculated. As 
M = {df + 9|)Aa; will not be used in the construction of Green function, we 
are only interested in the case where M is a homogeneous quadratic form in the 
component of Vx «• e. 

3 g2 



2 



i=i ^ 

So, the model equation (5.1) can be rewritten as: 

E-%-(i-/^-^H)^^;;^ -.^0, V, (5.2) 
j=i 3 

By a change of variables Xj = rnjrjj, equation (5.2) becomes Poisson equation 
in 'if{ri,uj) = 5'(x, w) i.e. 

A,* =(l-/3^(.))^^^^^ =<!>(,,.) (5.3) 
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where, 



T(r/) 



9 9 9 9 9 9 



+ 



= t(x) and $(7], = $(x, w) 



V "1 "2 "3 

Notice that the source $(77,0;) is symmetric with respect to elhpsoid r, i.e. 

$(r/,a;) =¥(t,w). 

Therefore, the solution 5* of the Poisson equation (5.3) is the potential field of 
a uniformly charged ellipsoid due to a charge density ^{t,uj). The potential 
field \1/ can be calculated with a classical approach using ellipsoidal coordinates. 
(See for example [24, 31] for the theory of potential problems in ellipsoidal 
coordinates.) 

For the solution of the Poisson equation (5.3) we recall following result from 
[31, Ch. 7, Scc.6]. 

Proposition 5.1 Let 



3 

^7 



2 — 1 and g{z) = {oLjh)^ + z 



and let Z{h, C.) be the largest algebraic root of f{z)g(z) = 0. Then the solution 
of the Poisson equation 

A'^YiC) = 47rx f 4 + 4 + 4) C e 

is given by 

/•OO 

y(C) = 2TTaia2a3 / x{h)I{hX)dh. 
Jo 

The integrand I{h, C.) is defined as 



1 



dz 



zdz Z >0 



'a \/giz) 

Hence, the solution of (5.3) can be given as 



Z <0 



27r6 



vl/(77,w) = — l-/3^(c.) — 



1 



-lK(uj)h 



-I{h, rf)dh 



or equivalently. 



*(x,w) 



1 - l3A{uj) 



g \/ —lK{uj)h 



I{h^ x)(i/l, m = 771177127713 



(5.4) 
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By a change of variable s = h ^z, we can write J(/i, x) as: 



/(/i,x) 



mh 



mh 



1 

s(/i,x) \/G{s) 
1 



ds h < T 



(5.5) 



VG{s) 



zds h > T 



with S{h,x) = h ^Z{h,x.) being the largest algebraic root of the equation 



where 



F{s)G{s) = 



Gis) = —g{h's)^n'^^,{V,{s)} 



(5.6) 



with Vj{s) = b'j + m'jS 
Remark 5.2 Note that, F{s) = corresponds to a set of confocal ellipsoids 



-1 2 



(5.7) 



such that t(x) = /i(0) i.e. S'(t) = 0. Moreover, S > if the ellipsoid h lies 
inside r and S < if the ellipsoid h lies outside r. 

5.2 Derivatives of the Potential field 

Now we compute the derivatives of the potential We note that /(/i,x) is 
constant with respect to x when h > t. So, 



dxk 



dS{h,K) 



1 



dxk ^G(5(/i,x)) 







for k ~ 1,2,3 and by consequence 
1 



or 



dxk 



dxk 



8pnm 



1 - PA{uj) 



oo g\/^K{i»;)/i 



h <T 



h> T 



dl{h, x) 



dxk 



dh 



1 - /3Aiuj) 



^ —lK{uj)h 



dS{h,x) 



dxk ^G(5(/i,x)) 



dh. (5.i 
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Now, we apply — — for I ^ 1, 2, 3 on (5.8) to obtain the second order derivatives 
oxi 

of*: 



-8/97r 



dxkXi 



1 - I3A{lo) 



dxi 



-lK.{ijj)h 



dS 1 



+ h-fiA(u) 



-1k:(cj)t 



dxk ^/G(S) 
dS{T) 1 



-.dh 



dxk ^G(5'(r)) j 
1 f d^S 1 dS dS G'{S) \ 
,/G{S) I dxkdxi ~2d^d^i G{S) J 



As F{S)G{S) = and G{s) is normally non-zero on S, therefore by differ- 
entiating F{S) = 0, we obtain [19, eq. (5.21)-(5.23)] 



-2xk 



dS 



dxk Vk{S)F'{S) 



(5.9) 



-AxkXi 



F"iS) , ml 



dxkxi Vk{S)Vi{S)[F'{S)f [F'iS) VkiS) 1^(5) J VkiS)F'{S) 



25. 



kl 



(5.10) 



where. 



F'{s) = y: 



2rnlx] 



G'{s) = G{s)J2 



(5.11) 

and prime represents a derivative with respect to variable s. 

Substituting the values from (5.9) and (5.10), the second order derivative of 
becomes 



4p7r 



dxkXi 



-XkXl 



(l - PA{^) 



aalafF'iO) 



\/ —lK(uj)h 



F'{S)^G{S) 
F"{S) ml 7nf I G'jS) 
Vk(S)ViiS)F'iS) { F'iS) + VkiS) ^ ViiS) + 2 G(5) 



2xkXi 



VkiS) 
(5.12) 



dh 



Remark 5.3 If for some i £ {1, 2,3}, — >■ one semi axis of the ellipsoid t 
tends to infinity but no singularity occurs. Therefore the results of this section 
are still valid in this case. 



14 



6 Elastodynamic Green Function 

In this section wc present the expressions for the elastodynamic Green functions 

for the media presented in section 3. Throughout this section Cr, — . with 

V P 



{1,2,--- ,6}. Werecah that Ki{uj) = _ 
6.1 Medium I 

Ah the eigenvectors of F are constants in this case i.e. T>i ~ e^, therefore 
Mi ~ 1 and =6^06^. If G is the Fourier transform of the viscoelastic Green 
function G with respect to variable t, then: 



^ 1 

G = ^$,(x,a.)e,®e, = — ^ 



c,+3 (l - pJiiu)'^ 



CiCiC^ceTi 



e,;(X>e,; 



(6.1) 



where 
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C4 


C3 








1 


Ml 


II 


Cl 


Cl 


C4 


1 


1 





M2 




C6 


C6 


C4 


* 


* 


* 


M3 




Cl 


Cl 


Cl 


1 


1 


1 


Ml 


III 


C6 


C6 


C4 


1 


1 





M2 




C4 


C4 


C4 


* 


* 


* 


Ma 



^2 ' 



Table 1: Values of bi and rrii for different media. Here * represents a value 
which is not used for reconstructing Green function. 
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6.2 Medium II 

According to section 4, the functions have following expressions: 

-lKi(ti;)ri(x) 



/ ^ \ p\/^'i2 (1^)1-2 (x) 

(x,a.) = (1 - /3,A{c.)) 



C4p7rT2(x) 

g\/^K3(;^)T3(x) 

4c|c4p7rr3(x) 



where 



n(x) = W-i + -| + -|, r2(x) = W-i + -! + -!, r3(x) = W-i + -| + -| 
V "^4 '^4 C3 y C;^ C4 y Cg Cg C4 

To calculate Green function, we use the expression 

G = $3! + Di (g) DiAff ^ ($1 - $3) + T>2® D2A/2~^ (*2 - $3) ■ 

Di = 63 and A/i = 1, yield 

Di (g) DiA/f ^ ($1 - $3) = ($1 - $3) 63 (g) 63 

To compute D2 (g D2A/2~^ ($2 - 'fa), suppose 

= A/2" ^$2 and *3 = A/2" ^$3 

and notice that mi = 7712 = 1 and 7713 = 0. Moreover for $2 and $3, 61 = 62- 
(See Table 1). Thus, we have 



Apt, 92*2 e e j ev^'=^(")-^ 1 1 , . 5 , 



RfcR; < 2 r Sol'^fci — 2RfcR, 







R-fcR-i < — 2 f T^^^ki — 2RfcR, 



rT3 




'0 





1 -/33yi(w)) "^fe^i [ CgC4T3 J C4/?^ 

where Rfc = — for k ~ 1, 2. See Appendix C for the derivation of this result. 
R 

By using the second derivatives of ^'2 and 4*3 and the expression 

2 

D2 ® D2A/2"^ ($2 - ^3) XI - *3) efc ® e, 

fc,Z=l 
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we finally arrive at 



gV'^K3(Lj)T3(x) , \ gV'-T«l(^)Tl(x) 



(1 - l32A{uj) 



4c|c4p7rr3(x) 

g%/^/i2(w)r2(x) 

4c^C4p7rT2(x) 



- 1 - kAiio] 



iclcsp-KTl (X 

gV^K3(tj)T3(x) 

4c|c4p7rr3(x) 



R® R 



(1 - 



(J - 2R® R)x 



d/i - 1 - psAiio) 



Or equivalcntly, 

G = $163 (g) 63 + $2R «> R + $3(1 - R «) R) 



1 



cf / h^2ih,uj)dh- Cq / h^3{h,uj)dh 
^0 



( J - 2R (g) R) 



Here J = I — 63 (g 63 and R = Riei + R2e2 
6.3 Medium III 

The solutions of the wave equation $i in this case are 

$2(x,w) = (1-/32^H) 



;(X,W) = (1-/33^H) 



4cfp7rri(x) 

g\/^K2((^)r2(x) 

4c^C4p7rT2(x) 
gV'^K;3 (1^)1-3 (x) 

Aclp-KTs (x) 



where 



n(x) = = -, T2(x) = + ^ + ^, T3(x) = - 

ci V ci y eg eg e| e4 

To calculate Green function, we once again use the expression 

G = $3! + Di ® DiMf ^ ($1 - $3) + D2 ® D2A/2"^ ($2 - $3) • 

Suppose = Mf ^$1 and ^'3 = Mf ^$3. Notice that nii = m2 = rns = 1 
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for Ml and 61 = &2 = &3 for $1 as well as $3 (see Table 1). Thus, 



4p7r 92^1 
dxkXi 



(1 - /3iA{uj) 



dxkXi 



(1 - f3,A{^) 




r 
















r'''3 




^— 1k3(cj)/i 


dh 


Jo 









See Appendix B for the derivation of this result. It yields 

r-l 



Di (8)DiMf ' ($1 - $3 
(1 - /3i^(w) 



1 

4p7r 



— 1ki(c»j)ti(x) 

cfri(x) 



+ 



(1 - PsA^) 



-lK3(i^)r3(x) 

C4T3(x) 



r (X) r 



1 - 



dh - [1- PaAiij] 



dh 



(I- 3f 0f) 



= - $3(x,a;)]f (8)f - — 



h^^i{h,uj)dh- / h^^3{h,uj)dh 



(I - 3r(g) r) 



where r = riei + r2e2 + ^363 with r ^ = — for alH = 1, 2, 3. 

r 

To compute, D2 D2Af2~^ ($2 - ^3), suppose *2 = M2^^2 and 'i'4 
M2^^3- By using formula (C.3) with mi = 7712 = 1 and 7713 = 0, we obtain: 



4p7r 92*2 
dxkXi 



(1 - /32-A(w) 



4p7r 9^*4 



(1 - /33-4:(o.) 




C4i?2 



C4i?2 



(4/ - 2RfcR,) 



(Jfc/ — 2Ra;R( 











Jo 












Jo 





dh 



dh 
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with Rfc ~ — and k,l G {1,2}. This allows us to write 
R 

D2 ® DaMj"^ ($2 - <^3) 



1 

Apir 



(1 - l32AiLj[ 



/^K2(cj)t2(x) 



1 - Mlj) 



ClT-2(x) 

(rIbi <S)ei - RiR2[ei ® 62 + 62® ei] + Rfe2 ® 62 



'^K3(w)r3(x) 

c1t3(x) 



4c4p7ri?2 



(1 - /32^(t^) 



/ — Ik. 2 (cj)/). 



- (1 - /33^(^)^ 

(1 - 2R^)ei (g) ei - 2RiR2[ei (g) 62 + 62 «) ei] + (1 - 2R?)e2 (g ea 









X 











= [$2 (x, w) - $3 (a;, w)] R^ (g R-^ 



1 



l"''2 rTa -| 

: / h<^2{h,uj)dh- cl / h^3{h,uj)dh ( J - 2R-L ® R 

where R''- = R2ei — Rie2 and J = I — 63 (g) 63. 
Finally, we arrive at 

G = $1? (g)f + $26.-'- (g R-*- + $3(1 - rig) f- R^ (g R-'-) 
^ \'^'S>i{h,uj)dh~ [ \^^z{h,Lo)dh 

n ./n 



1 



=1 / h^2{h,uj)dh ~ cl / h^3{h,uj)dh 
Jq Jo 



(I - 3rg) r) 

(J-2R^®R-^) 



6.4 Isotropic Medium 

When cge = C44. medium III becomes isotropic. In this case 

$2(x, = $3(x,i:^), ;Si=/32, ri(x) = — , and T2(x) = — = r3(x) 

Ci C4 

Thus, the Green function in an isotropic medium with independent elastic pa- 
rameters cii and C44 can be given in frequency domain as: 

G = $2! + Di (g DiMf ^($1 - $2) 

= $1? ®f + $2(1-? ®?) - A- / ^ h^^iih,Lu)dh~ / h'^^2{h,uj)dh (I-3f(gf) 

^ Jo Jo 

where <I>i and $2 are the same as in medium III. This expression of the Green 
function has already been reported in a previous work [18]. 
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A Decomposition of the Green function 

Consider the clastic equation satisfied by G: 

(r^(V.)G(x, t) + r (V.)^[G](x, t)) - P^^^ - S{t)Si^)I. (A.l) 
If G is given in the form 

3 

G = ^E,(V,)0, (A.2) 

i=i 

Then substituting (A.2) in (A.l) yields: 

sm^)i - (r(v.)G(x,i)+r(v.)^[G](x,t))-p^^=|^ 

= (^.^(V.)</'.+i}'(V.)^M)E,-(V.)E,,(V.)-p^E,(V.)^^^|^ 

By definition Ej(Vj;) is a projection operator which satisfies 

E,(V,)E/V,) ==%S,(V,) 
Consequently, we can have 

i.j — l i—1 

i—l ^ ^ 

3 

Moreover I = E— therefore 



E E.(v.) f (i?(v,)0, + m'^.)A[cf,,]) p^^^ - sm^)) = 



1=1 



Finally, remark that G we can express in the form (2.8) if the functions 
satisfy equation (2.10). 
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B Derivative of Potential: Case I 



If 61 = 62 = ^3 and mi = m2 = 7713, we have 
V,is) = V2is) = V3is)^bl + mls 



G{s) = {Vi{s)f and G"(s) = 6(5) 



3mf 



Viis) 



with r = y/xl+x^+x^. When F{S) = 0, we have 



1 



-1 



and 



-1 



2^ 
2 



TO^ 1 G'{S) \ _ 3 ?7if _ 3mlh 
F'{S) ^ Vk{S) ~^ Vi{S)^2 G{S) J ~ 2^1(5) ~ 2^ 

Substituting (B.l) and (B.2) in (5.12) wc finaUy arrive at: 



1 - f3A{uj)^ ^""kxi 



^ —1h.(uj)t 



(B.l) 



(B.2) 



dh 



(B.3) 



where r,- = — for aU j = 1,2,3. 
r 
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C Derivative of Potential: Case II 

If 61 =62, "^1 = TT^2 and m3 — 0, wc have 



Vi{s) = ¥2(3) = bl + mis andV3(s) = 6i 
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and F'(0) 



G(s) = 6^ (^1(5))^ and G"(.s) = G{s) 



(C.l) 



2m^ 



with R ~ \J^\ +2^2- For all fc e {1, 2}, we have 



1 



1 



-1 



V{S) 



and 



(C.2) 



F"{S) ml mf IG'jS) 
F'{S) Vk{S) Vi{S) 2 G{S) J Vi{S) 

Substituting (C.l) and (C.2) in (5.12) and simple calculations, we finally arrive 
at: 



= RfcR; < — 



1 



1 - /3^(w)) ^^fc^' 



br 



(Ski - 2RfeR0 







'0 






(C.3) 



where Rfe = ^ for fc = 1, 2. 
it 
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